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By carrying out Monte Carlo simulations based on the two-species atomic- 
scale kinetic growth model of GaAs(OOl) homoepitaxy and comparing the results 
with scanning tunneling microscope images, we show that initial growing islands 
undergo the structural transformation before adopting the proper /32(2 x 4) re- 
construction. 

In developing optoelectronic devices, GaAs(OOl) is often used as the basis sub- 
strate for fabrication HP . To study the atomic-scale growth kinetics of this surface, 
extensive use has been made of ah initio calculations with the energy-minimization 
procedurelJ. With this method, however, it is actually difficult to elucidate the 
true growth kinetics right at the growth temperatures, because growth proceeds in 
non-equilibrium conditions there. For example, by investigating the stable sites for 
Ga adatoms and assuming arsenic species to stick onto them, it is deduced in Ref.tf 
that the nucleation of islands is initiated by the sticking of Ga adatoms at the trench 
sites of the GaAs(001)-/?2(2 x 4) structure, which is depicted in Fig.l. If this growth 
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Fig. 1. GaAs(001)-/32(2 x 4) structure. The dark and bright disks denote, respectively, As and Ga 
atoms, whose radii decrease according to their depths. The shaded area denotes the unit cell. 

mechanism is true, a quite large portion of a GaAs(OOl) surface must reduce the 
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widths of the trenches from double to single before islands nucleate. However, no 
scanning tunneling microscope (STM) images obtained thus far seem to support it. 

Here we carry out atomic-scale Monte Carlo (MC) simulations as the alternative 
means to investigate nucleation and growth of islands on the GaAs(001)-/?2(2 x 4) 
surface. For this purpose, we developed the^ two-species model of homoepitaxial 
growth on the GaAs(001)-/32(2 x 4) structured. 




Fig. 2. Incorporation kinetics of As2 molecules, (a) Ga atoms are not properly arranged on a 
surface, so that the incorporation is prohibited. Only the desorption of As2 (R — » V) takes 
place, (b) Incorporation (R — > S) is enabled. The barrier we used for the R — * S process is 1.74 
eV, and the barrier difference between the two processes R — > S and R — > V is 0.65 eV. 

For the arsenic deposition process, we use As2 molecules, for which we simplify 
the incorporation process onto a surface structure in a mean field manner by intro- 
ducing the chemical reservoir, where deposited As2 molecules are once stored. This 
process is schematically drawn in Figs. 2 (a) and (b), where R, S, and V denote, re- 
spectively, the reservoir state, the surface-incorporated state, and the vacuum, which 
is far distant from a surface. 

As Figs. 2 (a) and (b) indicate, whether an As2 molecule can be incorporated 
onto a surface site depends on if surface Ga atoms are properly arranged there or 
not. When they are properly arranged, the incorporation channel of an As2 molecules 
becomes open, and the potential for As2 molecules changes from Fig. 2(a) to Fig. 2(b). 
As a consequence, the macroscopic growth rate is controlled by the deposition rate 
of Ga atoms. However, at atomic scales, arsenic species come into play in the growth 
kinetics, as we will see below. Note that this treatment of the chemisorption process 
improves the efficiency of computation, so that we can study the morphology of a 
growing surface in a realistic computation time. 

For the interactions between surface atoms, we take account of those between 
the same species in the same layers, as depicted for the Ga-Ga interactions in Fig. 3 
and for the As- As interactions in Fig. 4. 

Among them, the interactions in Figs. 4 (d) and (e), representing the Coulomb 
repulsions between two surface As atoms, play the central role in realizing the ,92(2 x 
4) structure u . We note that these interaction parameters were determined solely by 
the comparison between the calculation results and the STM images. 

By carrying out simulations with the kinetic MC algorithmic , we found that the 
nucleation of islands take place on top of the /?2(2 x 4) substrate without having 
the (2 x 4) structure, which transform into the (2 x 4) structure at a later stage of 
growth cP. This is because the excess charge on surface As atoms, which is the origin 
of the Coulomb repulsions in Figs. 4(d) and (e), is not significant on small islands. 
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Fig. 3. Ga-part of the interactions. The interaction energies are (in eV): (a)0.28, (b)0.25, (c)0.28, 
(d)0.32, and the back-bonding energy is 1.40 eV per atom. 




Fig. 4. As-part of the interactions. The interaction energies are (in eV): (a) 1.00, (b)0.63, (c)— 0.15, 
(d) — 1.21, (e)— 0.25. The back-bonding energy is 1.65 eV per atom. 

The processes of nucleation and the structural transformation are well exemplified 
by the series of the snapshots in Fig]B[ 

To further confirm the validity of our model quantitatively, we plotted in Fig. 6 
the island-size distributions obtained by carrying out the calculations on the 300 x 
300 lattice as well as on the 240 x 240 lattice and compared them with the STM 
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Fig. 5. Nucleation and growth of an island. 



images obtained at several different coverages. Here the lattice spacing in the model 
calculations corresponds to 4.0A on a real surface. The result of the comparison 
for the total number of islands is displayed in Fig. 6. In this comparison, the results 
of the calculations were obtained at T = 580° C, whereas the STM images were 
obtained after quenching the surface down to room temperature. In order to avoid 
the ambiguity caused by this difference in the conditions, we examined the effect 
of quenching, i.e., we performed the quenching simulations to find that single As 
dimers appearing at T = 580° C change into double ones by quenching, so that 
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Fig. 6. Total number densities of islands as functions of a coverage compared with STM measure- 
ments at the indicated coverages. 



the number of islands does not change during this process. We further found that 
other transient non-(2 x 4) structures tend to adopt the /32(2 x 4) structure during 
quenching. These results indicate that we can avoid the ambiguity and, hence, we 
can carry out a reliable comparison between the simulations and the STM images 
if we count the numbers of islands before, during, and after the splitting separately 
by classifying them into unsplit, partially split, and split islands, and discard the 
partially split ones, because they appear only as transient non-(2 x 4) structures. 
The results are seen in Fig. 7 to show a very good agreement between them. 

Since Ga atoms prefer to reside at trench sites, As2 molecules can stick onto them 
and the width of the trench does change locally from double to single, as suggested 
by the ab initio calculations BP. However, according to our calculations, these surface 
As dimers desorb from such sites within 1 ms at T = 580° C, so that actually these 
processes do not result in the growth of islands. Instead, when Ga atoms gather on 
top of the As-dimer row of the /32(2 x 4) structure and an As2 dimer sticks onto 
them, this is stable enough to be the nucleation site. Note that the complexity of 
this process is the origin of the significant delay in the nucleation time seen in Figs. 6 
and 7. 

Our results show that the nucleation of an island is controlled by Ga kinetics, 
whereas the structural transformation and, accordingly, the appearance of the /32(2 x 
4) structure on newly created islands is due to the splitting process seen in Figs. || 
(d-1), which are controlled principally by the Coulomb repulsions in Figs. 4 (d) and 
(e), and thus, by As kinetics. By further employing the rate equations, we confirmed 
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that these two are the only rate-limiting processes at the pre-coalescence regime of 
GaAs(OOl) homoepitaxyB*. 
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Fig. 7. Densities of split, unsplit, and partially split islands obtained by simulations on a 300 x 300 
lattice and STM observations. 

In summary, by carrying out atomic-scale growth simulations of GaAs(OOl)- 
/32(2 x 4) structure, we found that the nucleation of islands takes place on top of 
the (2 x 4) substrate and the islands grow via the structural transformation from 
non-(2 x 4) structure to the /32(2 x 4) structure. We further showed the validity 
of our model study quantitatively by comparing the results with those obtained by 
scanning tunneling microscopy observation. 
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